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ABSTRACT .  This  report  develops  a  general  spreading  loss  expression  and 
ray  tracing  procedure  for  use  in  sonar  detection  studies  particularly 
where  shallow  water  makes  bottom  topography  a  significant  factor.  Any 
analytically  describable  ocean  bottom  can  be  accommodated  by  these 
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FOREWORD 

Underwater  acoustics  is  a  field  of  great  complexity  about  which  little 
is  clearly  understood.  The  variability  of  the  ocean  medium  and  the  presence 
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Ocean-bottom  topography,  which  is  one  of  the  most  significant  environ¬ 
mental  factors  in  shallow  water  areas,  suid  its  effect  on  underwater  sound 
transmissions  are  examined  in  this  report.  The  analytical  techniques  pre¬ 
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INTRODUCTION 

Acoustic  detection  studies  often  rely  on  ray  tracing  procedures  which 
analytically  model  the  transmission  of  sound  in  an  underwater  environment. 
These  procedures  include  the  calculation  of  spreading  loss,  which  is  the 
resultant  decrease  in  sound  level  as  the  wave  front  spreads  over  a  gener¬ 
ally  expanding  area.  However,  most  applications  use  a  spreading  loss 
expression  which  omits  allowances  for  a  non-level  ocean  bottom  and  may 
lead  to  serious  inaccuracies  in  those  shallow-water  or  long-range  acoustic 
detection  studies  where  bottom-reflected  sound  is  of  major  significance. 

This  report  eliminates  such  trouble  by  developing  a  general  spreading  loss 
expression  and  ray  tracing  procedure  which  may  be  used  with  any  analytically 
describable  ocean  bottom. 

THE  COORDINATE  SYSTEM 

A  coordinate  system  using  orthogonal  axes  X  ,  Y  ,  Z  ,  and  angles  0 
and  0  as  shown  in  Fig.  1  is  used  throughout  the  derivation. 


FIG,  1.  The  Coordinate  System. 


BASIC  ASSUMPTIONS 

1.  The  oceeui  medium  is  made  up  of  a  series  of  horizontal  layers  of 
thickness  A  Z  ,  with  each  layer  containing  a  constant  velocity 
gradient  k  =  A  V/  a  Z,  where  A/  is  the  difference  in  sovind  veloc¬ 
ities  at  the  layer  boundaries. 

2.  Any  anedytically  describable  ocean  bottom  may  be  used. 

3.  The  ocean  surface  is  a  horizontal  plane. 
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Under  the  above  assumptions,  sound  will  travel  from  a  source  point 
Pq  to  a  point  Pjj  along  a  ray  path  which  is  contained  in  a  series  of 
vertical  planes,  with  each  vertical  plane  Tj  containing  the  ray  path 
between  the  and  the  j  +  1st  bottom  reflections.  T.  =  "^^-1 

and  only  if,  the  normal  to  the  bottom  at  the  ,1^  reflection  lies  on 
vertical  plane  ,  The  course  of  the  ray  path  on  each  vertical 

plane  will  be  determined  by  refraction  and  surface  reflection.  Refrac¬ 
tion  is  based  on  Snell's  Law:  /cos  6^  “  V-^j/cos  ,  where  , 

,  and  are  the  respective  sound  velocities  and  angles  of  incli¬ 
nation  of  the  ray  at  any  two  depths,  and  ,  on  that  -portion  of  the 
ray  between  successive  boundary  reflections. 

The  initial  ray  direction  can  be  uniquely  represented  by  two  para¬ 
meters,  Qq  and  0Q  ,  where  9^  is  the  angle  the  ray  initially  makes  with 
the  horizontal  and  0q  is  the  angle  between  the  vertical  Tq  and  ^Z 
planes.  Given  specific  environmental  parameters  and  a  ray  source  point 
Pq  ,  the  coordinates  (  X,  Y,  Z)  of  any  point  P  on  a  ray  are  functions 
of  initial  ray  direction  ,  0  )  and  some  third  parameter  such  as 
travel  time,  ray  path  length,  or  horizontal  distance  covered,  which  fixes 
the  location  of  the  point  on  the  ray. 

THE  SPREADING  LOSS  TERM 

The  spreading  loss  term  (SL  )  is  the  ratio  of  acoustic  intensity 
at  point  Pjj  to  the  acoustic  intensity  at  an  index  point  a  unit 
path  length  from  source  Pq  .  All  acoustic  energy  is  contained  within  a 
ray  bundle  as  defined  in  the  next  paragraph,  and  a  sound  wave  front  at 
creates  a  ray  bundle  cross-sectional  area  If  energy  is  assumed 

evenly  distributed  over  A^^  ,  it  follows  that  I^  is  inversely  proportional 
to  A  .  Therefore,  if  A^  is  a  similar  cross-sectional  area  at  the 
index  point, 


(3^)n  =  ^n/^o 


(1) 


RAY  BUNDLE  CONCEPT 

Given  a  basic  ray  which  leaves  P^^  with  initial  direction  (  6o»^o 
a  ray  bundle  is  that  volume  which  is  bounded  by  the  four  ray  paths 
leaving  Pq  with  initial  directions: 

1.  («„ 

-‘‘V 
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where  AQ  and  A  0  are  infinitesimal, 
o  o  , 

Bo\inding  rays  1,  2,  3,  and  2+  can  be  generated  by  considering  two 
other  rays,  5  and  6,  whose  respective  initial  directions  are  +  A 
0Q  )  and  (Oq  -  ®o»  ^o^  shown  in  Fig.  2.  These  two  rays  and  the 

basic  ray  initially  lie  on  the  same  vertical  plane,  Tq  .  By  rotating 
Tq  through  the  incremental  angle  A0q  about  the  vertical  axis  passing 
through  Pq,  rays  5  and  6  form  rays  1  and  3,  respectively,  which  initially 
lie  on  a  common  vertical  plane,  .  By  rotating  Tq  through  the  angle 
-/\0Q,  rays  5  and  6  form  rays  2  and  4,  respectively,  which  initially  lie 
on  a  common  vertical  plane,  T^’. 


P. 


FIG.  2.  The  Ray  Bundle. 

CROSS-SECTIONAL  AREA  OF  A  RAY  BUNDLE 

The  ray  bundle  cross-sectional  area  A^  is  defined  as  the  area  of 
the  soimd  wave  front  contained  in  the  bundle  at  point  Pn  on  the  basic 
ray.  A  popular  approximation  of  A^  is  obtained  from  the  product  of 
equations  3B-38  and  3B-40  in  Ref.  (1) ,  where  it  is  assumed  that  the  basic 
ray  and  the  bounding  rays  are  contained  intirely  in  their  respective 
initial  vertical  planes.  However,  since  this  assumption  is  not  valid 
after  the  occurrence  of  a  reflection  from  a  sloping  bottom,  a  more 
general  approximation  of  A^  must  be  derived  before  any  non-level  ocean 
bottom  is  considered. 

Let  tn  be  the  time  required  for  a  sound  wave  to  travel  from  source 
Pq  to  point  Pjj  .  Ajj  is  the  area  of  that  surface  which  is  described  by 
the  trace  of  point  Pjj  as  the  initial  ray  direction  varies  from  (0^  ,0q) 
to(flQ  ±  ^  ©Q,  0o  *  ^  0o  "^n  remaining  constant.  Since 

and  a0q  are  very  small,  k-^  is  a  function  of  Pjj  ,  aPj^/ aSg  ,  and 
3  0q,  where  atj^  /  ^  ©0=  0  and  ®  tj^/  ^  0o  =  0. 

The  first  basic  assumption  indicates  that  the  ocean  medium  has  a 
velocity  structure  in  which  the  velocity  gradient  is  a  piecewise  con¬ 
stant  function  of  depth.  The  coordinates  of  point  P^  are  (X^  ,Yjj  ,Zjj  ) 
where  Zj^,  the  depth  of  P^j,  may  be  the  depth  of  a  velocity  gradient 


3 


NAV^'fEPS  REPORT  7799 


discontinuity.  Therefore,  it  is  conyenient  to  approximate  from  P  , 

aPn/30o  and  with  aZn/^  ®o  =0  and  3Z^  1 3  0^  =  0  . 

Consider  the  quadrilateral  Qj^  with  vertices  P^^j^  ,  ^2-^  «  P3n»  and 

Plijj  which  are  the  respective  points  of  intersection  of  bounding  rays 
1  through  4  with  horizontal  plane  passing  through  point  P^  .  The 
area  of  is  obviously  the  horizontal  cross-sectional  area  of  the  ray 
bundle  at  P^  .  The  desired  cross-sectional  area  is  approximately 
equal  to  the  area  of  the  perpendicular  projection  of  onto  H*  ,  the 
basic  ray's  normal  plane  at  Pj^.  If  o(^  is  the  angle  between  planes 

lln  'in'i  ''-Jl  . 

A^  =  cos  cL  ^  (area  of  Q^j).  (2) 

An  example  of  the  perspective  of  planes  11^  and  H*  is  given  in  Fig.  3. 
The  ray  direction  at  point  P^j  is  (©^  ,  0jjj )  »  where  is  the  ray's 
angle  of  inclination  and  0^^  is  the  angle  between  the  vertical  XZ  plane 
and  the  vertical  plane  containing  the  ray  path  after  the  mth  bottom 
reflection. 


It  can  be  seen  from  Fig.  3  that  <=<^^=90°-  0^.  Substituting  this  into 
equation  (2), 

A.,  =  sin  ©„  (area  of  C-L).  (3) 

The  basic  ray,  defined  with  initial  direction  (©^  ,  0^),  intersects 
horizontal  plane  at  point  P^,  which  is  at  depth  Zj,  .  Rays  5  and  6, 
previously  defined  with  respective  initial  directions  (  ©q  +  ^  ©^  >0o  ) 
and  (  ©^  -  A  ©^  ,  0p  ),  will  Intersect  Hjj  at  points  P^j^  and  Pg^^,  where 
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P5n  =  Pn  +  <^Oo)  ^  ©o 

with  ^2^/  3©^  =  0. 


P6n  =  Pn  -  (»Pn/®eo) 

(1^) 


Bounding  rays  1  through  ^4,  defined  with  respective  initial  directions 
(©  +Ae^  ,  5^0  ^  ).  (®o  »  5^0  -  ).  (Qo  -  -^^o  + 

),  and  (Qq  -  ^©0  »  5^0  “  intersect  Hn  at  points  Pin 

through  Pl^n  *  the  vertices  of  quadrilateral  Qn  •  Therefore,  using 
expression  (I4),  the  vertices  of  are 

Pin  =  P5n  +  (9P5n/®«*o)^?'o  P2n  =  P5n  -  (aP5n/a0o)  "" 

Pan  =  P6n  +  (^P6n/3  0o)^«>o  P^n  =  P6n  "  (3P6n/^«>o)  K  ^5) 

with  ^  ^n^^  0o  “ 


Pljj  through  Pi|fj  can  also  be  expressed  in  terms  of  two  other  rays,  7 
and  8,  whose  respective  initial  directions  are  (©q  ,  0o  ‘*'^0o  ^ 

(®o  »  0o  ~  ^  0©^  •  7  and  8  will  intersect  Hjj  at  points  P^^^  and 

P8n  »  '"here 

P7n  =  Pn  ^  (^Pn/^0o)"^0o  =  Pn  >  ^  ®Pn/ 

with  az^/a0o  =  0. 

Using  expression  (6),  the  vertices  of  are 

■’in  =  ^Tn  *  Pan  '  -"Sn  *  OP6n/3ec) 

P3n  =  Pin  -  <®P7n/a»o)*»o  Pl.„  =  '  <®  P8n/3«o!^8„  '' 

with  aZn/ae^  =  0. 


Expressions  (5)  and  (7)  indicate  that  P^jj  ,  Pg^  ,  Pvjj  ,  and  P3JJ 
are  midpoints  of  the  sides  of  Qjj  ,  and  expressions  (4j  and  (6)  indicate 
that  Pn  is  the  point  of  intersection  of  lines  Joining  opposite  midpoints. 
Figure  4  is  a  diagram  of  with  the  related  points  described  in  ex¬ 
pressions  (4)  through  (7). 


It  is  proved  in  the  Appendix  that  the  area  of  any  quadrilateral 
equals  eight  times  the  area  of  any  triangle  whose  vertices  are  the  mid¬ 
points  of  two  adjacent  sides  of  the  quadrilateral  and  the  poin^  of 
intersection  of  the  two  lines  Joining  opposite  midpoints.  If  is  the 
area  of  triangle  PnP5nP7n 

Area  of 


5 


NWWEPS  REPORT  7799 


Combining  equations  (3)  and  (8), 

^n  =  si"  ®n 

The  coordinates  of  point  are  (X^j  ,  Zj^).  Therefore, 

expressions  (4)  and  (6), 

Coordinates  of  P^j^  =  +AX5JJ,  Y^  +AY5JJ,  Z^) 

Coordinates  of  P^^  =  (X^^  ^n  ^^^Tn*  '^n) 

where 

AX5n  =  {9XJ3QJ£.Q^  AY5„  =  (aYn/^eo)^©^ 


The  relative  positions  of  points  Pn  ,  P^jj  ,  and  Py^j  are  shown 
5»  where  the  coordinate  system  origin  is  at  P^. 


(9) 

rom 

(10) 


(11) 


in  Fig. 
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Y 


FIG.  5.  Triangle  PnP5nP7n. 


From 


Fig.  5,  the  area  of  triangle  PnP5nP7n 

'i„  i^5niY5n  ^5n 

”  *^5n“^Tn - ^ - ^ 

AX7nAY3n  - 


2 


Combining  equations  (9)  and  (12), 


^X7„)(AY7n  -  AY3„)1 

2  j  (12) 


=  It  sin  (AXyjjAY^n  (13) 

Substituting  expression  (11)  into  equation  (13),  the  desired  expression 
for  the  ray  bundle  cross-sectional  area  at  point  Pn  is 

/3Xn  3Y  SX  9Y  \ 

\  Bin  8„  -3^  S^)  (IM 

There  is  a  certain  advantage  in  retaining  the  algebraic  sign  of  A^ 
even  though  it  has  no  significance  as  applied  in  equation  (l).  If  A^j 
and  Ajj+i  have  opposite  signs,  the  ray  bundle  has  passed  through  a 
focussing  point  somewhere  between  points  and  P^+i  • 

RAY  BUNDLE  CROSS-SECTIONAL  AREA  AT  THE  INDEX  POINT 

Ap  is  the  sound  wave  front  area  in  the  ray  bundle  at  the  index 
point  unit  path  length  from  the  ray  source.  Assuming  that  the  ray  path 
is  a  straight  line  over  the  first  unit  length,  Aq  is  approximately  the 
area  of  that  surface  which  is  generated  by  the  movement  of  the  index 
point  as  the  initial  ray  direction  varies  from  (q^  ,0q  )  to  (Oq  +a©  , 

0o  ±  ^  )»  S’®  described  in  Fig.  6. 
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FIG.  6.  Area  A  . 

o 

As  ©o  is  varied  from  -aQq  to  0  +  /^  ,  the  index  point  de¬ 
scribes  an  arc  of  length  2^0^  ,  lying  on  vertical  plane  T  .  As  0^ 

is  varied  from  0q  0q  to  0^  +-A  0q  ,  the  index  point  describes  an 
arc  of  length  2a 0^  cos  0^  ,  lying  on  a  horizontal  plane.  Therefore, 
the  ray  bundle  cross-sectional  area  at  the  index  point  is 

Aq  =  4a0qA0q  cos  ©0  (15) 


AN  OCEAN  BOTTOM  BOUNDARY  CONDITION 

Ray  bundle  cross-sectional  area  A^j  was  derived  from  points  , 

,  and  P7n  ,  the  respective  intersections  of  the  basic  ray  and  rays 
5  and  7  with  horizontal  plane  Hjj  .  However,  when  P  is  the  intersection 
of  the  basic  ray  and  a  non-level  ocean  bottom,  points  and  are 
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not  uniquely  defined  since  rays  5  and  7  will  Intersect  plane  Hjj  twice. 
This  difficulty  is  circumvented  by  assuming  that  points  and 

apply  to  rays  5  and  7  before  reflection  and  by  introducing  two  new  points, 
P^jj  and  P|j^  ,  which  apply  to  the  rays  after  reflection.  Two-dimensional 
examples  of  the  situation  are  shown  in  Fig.  7  and  8. 


Z  2 

FIG.  7.  Up-Slope  Reflection.  FIG.  8.  Down-Slope  Reflection, 


The  coordinates  of  point  P^  are  (  X^,  Yj,,  Zjj).  When  P^^  is  a 
bottom  reflection  point,  the  before-  and  after-reflection  rates  of 
change  of  Xjj  and  shall  be  respectively  represented  by  unprimed  and 
primed  partial  derivatives.  Expressions  (lO)  through  (lU)  will  therefore 
determine  the  before-reflection  ray  bundle  cross-sectional  area,  with  the 
after-reflection  area  similarly  found  from  points  P^,  ,  and  P|jj 

where 


Coordinates  of  P^^^  =  (Xj,  +AX^n»  ^^^n*  ^n^ 

Coordinates  of  P|j^  =  (Xjj  +.AX|jj,  Yj^  + 


with 

AYj„  =  (aYy50o)^0o 


(17) 


(18) 


It  remains  to  relate  the  above  after-reflection  terms  to  the  before¬ 
reflection  terms  in  expressions  (lO)  and  (ll),  thereby  obtaining  the 
primed  partial  derivatives  of  X^^  and  Y^  which  may  be  used  in  equation 
(l4)  to  calculate  the  after-reflection  cross-sectional  area. 

For  the  moment,  let  Pjj  be  the  basic  ray's  mth  bottom  reflection 
point.  Let  the  ocean  bottom  tangent  plane  at  this  reflection  point  be 
Ej^  ,  defined  by  and  two  angles  of  inclination,  and  which 
are  the  angles  plane  Ejjj  makes  with  the  horizontal  in  the  respective 
XZ  and  YZ  planes.  Let  the  bottom  reflection  entrant  and  emergent 
basic  ray  directions  be  (O^  ,  0in_i  )  and  (0^  ,  0ui) ,  where  0  is  the  ray's 
angle  of  inclination  and  0  is  the  angle  between  the  vertical  XZ  plane 
and  the  vertical  plane  containing  the  ray  path.  The  respective  initial 
directions  of  the  basic  ray  and  rays  5  and  7  were  previously  defined  as 

(Qo  *^0  »  **0  *  ^o  ^o 
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A0q  are  infinitesimal,  adequate  approximations  of  desired  points 
and  are  obtained  by  assuming  that  rays  5  and  7  travel  in  straignt 

lines  with  basic  ray  entrant  direction  (  )  from  respective 

points  Pcjj  and  P^^  to  plane  ^  ,  reflect,  and  travel  in  straight  lines 
with  basic  ray  emergent  direction  (  )  from  to  points  P^^^  and 

.  An  example  of  the  relationship  of  point  P^j^  to  Pjj  and  P^^  is 
shown  in  Fig.  9*  where  plane  Ejj^  is  described  by  point  P^j  ,  a  negative 
and  a  positive  cr^.  (The  relationship  of  P^^  to  P^  and  Py„  will 
not  be  discussed  further,  since  it  is  similar  to  that  for  P^^j  .) 


FIG.  9.  Relative  Positions  of  Points  ,  and  p^^. 

From  Fig.  9,  the  vertical  distance  between  noint  P5n  and  plane 
is 


Z  =  tan  tan  CT-^. 

A  top  view  of  Fig.  9  is  shown  in  Fig.  10  where,  noting  that  and 

AX'  are  negative. 


AX'  =  AX^  +  X  +  X' 

5n  5n 

AYl  =  AYc  +7+1' 

5n  pn 


(20) 
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A  closeup  of  ray  5  as  it  travels  fron  point  P5n  to  P'  is  shown  in 
Fig.  11,  where  it  can  be  seen  that 


X  =  M  cos  =  Z>  cos 

tan  0, 


n 


X'=  M'  cos  0„  =  • 

-  Z' 

cos 

lu 

tan  0^ 

Y  =  M  sin  0_  1  = 

Z’ 

sin 

®^in-l 

lU^JL 

tan  0„ 

(21) 


Y*=  M'  sin  0JJ,  = 


Z’  sin  0 


tan  0* 
n 


m 


Substituting  expression  (21)  into  (20), 

^X'  =  AXc  +  Z*(cos  0„  ,  tan  0'  -  cos  0„  tan  0^) 

tan  0  tan  O' 
n  n 

^Y^n  =  Z'(sin  0n,.i  tan  -  sin  0,,  tan  0^) 

tan  0  tan  O' 
n  n 

Further  inspection  of  Fig.  11  yields 
Z"'  =  Z  +  X  tan  '/'nj  +  Y  tan<r^ 

=  Z  +  X(tan  '('jj  +  tan  0j^_i  tan  (T"^) 


(22) 


(23) 
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Combining  equation  (23)  and  the  first  equation  in  expression  (2l) , 

Z'  =  Z  +  Z*  cos  (tan  'I' +  tan  tan 

tan  Qn 


=  Z  +  Z*  (cos  0^  ^  tan  +  sin  tan 


-  Z  tan  Qn 

“  tan  -  cos  0jjj_i  tan  <1*  ^  -  sin  0^.1  tan  er'^ 
Substituting  the  above  into  expression  (22), 
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s. 


-  ^>^50  *  <=ml  '^''sn  %  *  ^hn 

■*"  ^in2  m  ‘‘’  *^m) 

Substituting  pertinent  equations  from  expressions  (11)  and  (l8)  into 
the  above,  the  desired  ocean  bottom  boundary  condition  compensation  for 
ray  5  is 


- 

ax„  + 

^ml 

/  aXjj  tan 

^  m  + 
^  HI 

1 

o 

CD 

<D 

O 

^Oo 

+ 

^m2 

/a 3^  tan 

(J;  „  + 
^  HI 

^®o  " 

5^0 

^Qo 

(26) 


Expression  (26)  displays  a  relationship  of  terms  obtained  from  points 
Pn*  P5ni  ^^n*  Since  the  discussion  on  page  9  indicates  that  rays 
5  and  7  are  assumed  to  follow  parallel  courses  (both  v/ith  basic  ray 
direction)  in  the  space  immediately  adjacent  to  the  basic  ray’s  bottom 
reflection,  expressions  (lO),  (ll),  (17),  and  (l8)  are  used  to  acquire 
the  similar  relationship  of  corresponding  terms  from  points  Pn,  P7n» 
and  P|jj: 
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a^^n  ^  C^l 

/aXj^  tan 

'I'm 

+  3Y 

n 

300 

w 

3  00 

^^n  ^  C„2 

( 3'L^  tan 

<!>  m 

-t  a^n 

a5^o  ( 

3  0g 

tan 


n) 


tan  <7-_ 


(27) 


In  order  to  uniquely  identify  tottom  reflection  points,  let  n  =  Bm 
at  the  mth  bottom  reflection.  Then  Pg^j^  is  the  point  of  reflection,  with 
(©Bjjj,  and  the  bottom  reflection  entrant  and  emergent 

basic  ray  directions.  The  ocean  bottom  tangent  plane  at  this  reflection 
point  is  Ejjj,  defined  by  point  Pg  and  two  angles  of  inclination,  i|;  and 
a~  ,  which  are  the  angles  plane  %  makes  with  the  horizontal  in  the 
respective  XZ  and  YZ  planes.  The  coordinates  of  point  Pqjj,  are  (Xgjj^,  Yjjj^, 
ZBni).  Given  5  Qg ,  aYgm/ ,  and  aXBm/^0o.  and  aYsm/S^o 

immediately  before  reflection,  expressions  (25)  through  (27)  yield  the 
following  corresponding  terms  immediately  after  reflection; 


+  C^i  /  aXg^  tan  tan 


3©o 

a©o  ^ 

aeo 

aYSm 

^^Bm  Cm2  1 

(  9X3^  tan 

'1'  m  ^^Bm 

^®o 

a©0  1 

a 

where 


air 

OY* 

Bm 

TJT 

^ml  = 


a^Bm  +  Cmi  I  ^XBm  tan  <1^  +  SYgm  tan  cT-nj 

Wo . 


a  0c 


3  00 


S-^Bm  ^  Cm2 


^%m  '‘'m  *  tan  ff-j, 


5  0o 


P  0c 


(28) 


cos  <^m-l  tan  -  cos  0^,  tan  O3  _ 

tan  0;^^  (tan  ©Bm  "  cos  0^_i  tan  ij; -  sin  0j„_i  tan  cTj^) 


Cj„2  =  s^n  0n-.-l  tan  -  sin  0^^  tan  ©Bm _ _ 

tan  QJm  (tan  ©g”  -  cos  tan  i);  m  -  sin  tan  <7“^) 


Let  (SL)Bni  ^•nd  (SL)Jjjj  be  the  respective  spreading  loss  terms 
immediately  before  and  after  the  mth  bottom  reflection.  Then,  from 
equation  (I6), 
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cos  ©Q 

/  3XBm 

^^Bm  - 

S^Bm 

^^Bm 

sin  ©Bm 

U©o 

300 

3  00 

cos  ©Q 

•3Yim  - 

3X^m 

sin  ©^m 

U^o 

3  00 

aOo  . 

-U 


-11 


(29) 


where  the  -primed  partial  derivatives  are  obtained  from  expression  (28). 
CALCULATION  OF  SPECULARLY  REFI.ECTED  RAY  DIRECTION 


The  mth  bottom  reflection  entrant  and  emergent  ray  directions  are 
(©Bm*  ^m-1^  (®fim»  ^^m) »  where  ©  is  the  ray's  angle  of  inclination  and 
0  is  the  angle  between  the  vertical  XZ  plane  and  the  vertical  plane. con¬ 
taining  the  ray  path.  It  is  desired  to  express  and  0j„  in  terms  of 
©Bm  and  0ni_i. 

Let  Vjjj  and  be  the  respective  mth  bottom  reflection  entrant  and 
emergent  ray  tangent  unit  vectors.  Then,  using  Fig.  1,  the  direction 
cosines  of  and  are 

Vjjj  =  (cos  ©Bn  cos  0nj_i,  cos  ©Bm  sin  0ni.i,  sin  ©Bm) 

(30) 

=  (cos  ©§m  cos  0n,  cos  sin  sin 

Let  Nj^  be  a  unit  vector  which  is  normal  to  the  ocean  bottom  at  the 
reflection  point.  Njjj  is  also  normal  to  Ej^^,  the  bottom  tangent  plane 
defined  in  the  preceding  section  by  the  point  of  reflection  and  two 
angles  of  inclination,  and  Therefore,  the  direction  cosines 

of  Nm  can  be  expressed  in  terms  of  m  and  Figure  12  is  an  example 

of  bottom  normal  N^j,  with  bottom  tangent  plane  Em  described  by  a  negative 
and  a  positive  Noting  from  Fig.  12  that  the  direction  cosine 

ratios  of  Nm  are  (tan  ij^m*  "tan  (7"m*  “^)»  i'*'  follows  that  the  direction 
cosines  of  Nm  are 


Nm  =  /^ tan  it>  m»  tan  (T m,  -1  \ 

\yr  yr  yr) 

where 

C=  tan^  '^m  ‘^m  ^ 


(31) 
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Specular  reflection  imposes  two  conditions  on  the  entrant  and 
emergent  ray  tangent  vectors,  and  Vj?j,  and  the  bottom  normal, 

These  conditions  are 

=  aVj^  +  bNni  (V^i,  VjJi,  and  Nm  are  coplanar.)  (32) 

V^*Nm  =  -Vjn'Nm  (Hni  bisects  angle  between  -Vjjj  and  V^.)  (33) 


Since  the  above  vectors  were  defined  with  unit  lengths, 


V„*V_  =  V'*V'  =  =1 

mm  mm  mm 


(31*) 


It  remains  to  determine  from  equations  (32),  (33),  and  (3l*)  a  solution 
for  VjJj  in  terras  of  and  N^j,  which  may  be  used  with  direction  cosine 
expressions  (30)  and  (31)  to  solve  for  ^Bm  and  From  equation  (32), 


V.?i-(V: 


m 


VA)  = 


aVi„-(V„i  +V^)  +  bl>Ij„-(V^  + 


VA-Vm  +  V;-V^  =  a(Vj,-V^  +  V^.V^)  +  b(N^,.V^  + 
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Substituting  equations  (33)  and  (3^)  into  the  above, 
+  1  =  a(l  + 

a  =  1 

Equation  (32)  now  becomes 


v;.=  Vj,  +  bN„ 


From  equations  (33),  (3^),  and  the  above, 


b  =  -2Vni.Nm 


Substituting  the  above  into  equation  (35), 
is 


(35) 


the  desired  solution  for 


V;  =  -  2(V„.N^)N^  (36) 

Turning  to  the  application  of  equation  (36),  note  that  the  dot 
product  of  two  unit  vectors  is  equal  to  the  sum  of  the  products  of  their 
corresponding  direction  cosines.  Therefore,  from  expressions  (30)  and 
(31), 


Vm*Nm 


V  'N 


cos  Qbjjj  cos  0j„_i  tan  +  cos  ©Bm  sin  0^-1  tancT^jp  - 
sin  Qsm  ^  ^ 

Yc 

ys"  cos  Ogni  ^ cos  0,n,i  tan  »  ^  sin  0p,i  tan  CT-^  -  tan  Qbi^ 
C^3/T'  cos  (37) 


where 


Cin3  '  cos  0,n_i  tan  +  sin  tan  -  tan  Ogm 

tan^  tan^^m  ^ 

Using  expressions  (30),  (31),  and  (37)  with  (36),  the  direction  cosines 
of  are 
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COG 


cos  0j„  =  cos  ©Bjn  cos  -2 (0^13  cos  ©3^1)  tan 

Tr“ 


=  cos  ©B^  (cos  0JJ1.1  -  2C„3  tan  1!)^^^)  (38) 

cos  ©^j^  sin  0jn  =  cos  ©Bm  (sin  0m_i  -  20^3  tan  <r-^)  (38a) 

sin  ©;^ni  =  sin  ©Bn,  +  2Cn,3  cos  ©b^ 


The  above  equations  immediately  lend  themselves  to  the  following  rela¬ 
tionships  between  the  mth  bottom  reflection  entrant  and  emergent  ray 
directions,  (©Bm»  ^m-l^  (®5m»  ^m)  •  • 

©^n  =  sin  ©Bm  +  2Cm3  cos  ©Bm  (39) 

tan  0m  =  ^in-1  "  ^Cm3  . 

CCS  0„_3L  -  2C^3  tan  ij;  ^ 


=  cos  0p'_j,  tan  +  sin  0m_i  tan  -  tan  ©gm 
tan^  'I' m  tan^  cT'm  1 

EVALUATION  OF  THE  GENERAL  EXPRESSION  FOR  THE  SPREADING  LOSS  TERM 


(>♦1) 


Restating  equation  (I6),  the  general  expression  for  the  spreading 


loss  term  at  any  point  on  the  ray  path  is 


(SL)n  = 


-11 


cos  ©Q  /sXjj  __  <^^n  ^^n 


sin  ©n  (aSo  500  ^0o 


{h2) 


It  remains  to  derive  the  stated  partial  derivatives  of  Xj^  and  Y^. 


Point  Pj^  has  coordinates  (Xjj,  Yjj,  which  are  functions  of  initial 
ray  direction,  (Qq,  0o) »  ^^d  some  third  parameter  fixing  the  location  of 
the  point  on  the  ray.  Let  this  third  parameter  be  S^j,  the  total  horizontal 
distance  covered  by  the  ray  from  source  Pq  to  point  P^.  It  was  noted 
after  the  basic  assumptions  that  a  sound  ray  path  will  be  contained  in 
a  series  of  vertical  planes,  with  each  plane  containing  the  path  between 
successive  bottom  reflections.  A  top  view  of  a  ray  path  with  m  bottom 
reflections  is  shown  in  Fig.  13,  where  AS^  is  the  horizontal  distance 
covered  as  the  ray  travels  from  the  mth  bottom  reflection  to  any  point 
Pji  before  another  bottom  reflection,  and  0  is  the  angle  between  the 
vertical  XZ  plane  and  the  vertical  plane  containing  the  ray  path. 
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Horizontal  Distance. 


X 


4.) 


The  coordinates  of  the  ray's  mth  bottom  reflection  point  are  (Xgj^,  Ygj^, 
Zqjjj).  Therefore,  from  Fig.  13, 


^n  =  ^Bm  ^os  (2^, 

^n  =  Ygm  +  ASjj  sin  0^ 


(43) 


The  boundary  condition  discussion  preceding  expression  (28)  used  un¬ 
primed  and  primed  partial  derivatives  to  respectively  denote  before- 
and  after-reflection  rates  of  change.  Bearing  this  in  mind,  expression 
(43)  yields 


:  5X5m- 

A 

Sn 

sin 

o 

CI> 

300 

3Xn  - 

A 

Sn 

cos 

aOo 

9^n  _ 

,  - 

A 

Sn 

sin 

■ 

sio 

aYn 

9(Zlo 

A 

Sn 

cos 

300  .  5  00 

K  K 

900  ^©0 

^m  +  cos  0„  3A®n 


(44) 


where  the  primed  partial  derivatives  are  found  from  3Xgj^/  30o, 
9YBm/9«o»  9Xgjjj/a0o,  and  aYfim/ 90o.  using  expression  (28). 


A  general  ray  tracing  procedure  now  begins  to  take  shape.  Assuming 
that  expressions  for  ASjj,  9^3^/300,  90^^ 

5 0m/ 9 00  are  available,  the  spreading  loss  term  at  each  desired  point 
on  a  ray  is  obtained  from  expressions  (42)  and  (44)  where,  before  the 
first  bottom  reflection. 
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m  =  0 

300  =  0 


a 


.  0 


J  30  =1  Bm  =  0 

ni.  o 

3X^300  =  0  3X^300  =  0 


m 


a0„/30  =0  (45) 

m  o 

(since  0q  and  0^  are  independent  initial  parameters.) 

When  the  ray  reaches  the  mth  bottom  reflection  point,  0^^  =  Sgj^, 
3X„/a0o=  3X3^30^  3V3eo=  3W30O 

(46) 

3V30o=  3V30o=  3W30O 

While  at  the  mth  bottom  reflection,  expression  (44)  is  up-dated  by 
finding: 

1.  0j^  from  equation  (4o). 

2.  30  /30r,  and  3  0_/30^.  (method  not  yet  shown) 

j’fj  w  m  o 

3 .  3  XAm/  3 0o ,  SYiJ  30^ ,  3Xj^^/  3  0^ ,  and  ^Y^J  ^  0^  from 
expression  {28). 

V/hen  the  ray  leaves  the  mth  bottom  reflection,  0^  =  Qgjjj  for  an  instant, 
where  0^^^^  is  found  from  equation  (39)  •  After  reflection,  spreading 
loss  terms  at  desired  points  on  the  ray  are  obtained  from  expression 
(42)  and  up-dated  expression  (44). 

Expression  (44)  requires  the  evaluation  of  ^Sj^,  previously  defined 
as  the  horizontal  distance  between  the  mth  bottom  reflection  and  point 
Pjj.  Let  be  the  horizontal  distance  covered  as  the  ray  travels 

between  two  successive  depths,  and  .  Since  and  Z^  are  the 

respective  depths  of  the  mth  bottom  reflection  and  point  Pjj, 


n 


i=l  +  Bm 


(47) 


In  order  to  derive  AS^_j^,  recall  that  the  ray  path  is  refracted  accord¬ 
ing  to  Snell's  Law: 


Vg^/cos  0g  =  V^/cos  O3  between  successive  boundary  reflections 
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where  Vg^,  V-jj,  0^,  and  are  the  respective  sound  velocities  and  angles 
of  inclination  of  the  ray  at  any  two  depths,  and  Z-j^.  Since  the  ocean 
surface  is  assumed  to  be  a  horizontal  plane,  the  ray's  entrant  and 
emergent  angles  of  inclination  at  a  surface  reflection  will  be  of  equal 
magnitude.  Therefore,  the  preceding  equation  can  be  extended  to  read 

Va,/cos  ©g  =  V-^j/cos  ©^  between  successive  bottom  reflections,  (U8) 

Equation  (^+8)  invites  an  expression  of  velocity  as  a  function  of  ©  and 
Vy,  where  is  the  reversal  velocity  —  that  velocity  at  which  the  ray 
reverses  vertical  direction  because  of  refraction. 

V  =  Vp  cos  ©  between  successive  bottom  reflections  (^9) 


where 


Vr  =  Vg/cos  ©g  =  V^/cos  Q^. 

From  the  first  basic  assumption,  the  ocean  velocity  structure  con¬ 
sists  of  a  series  of  horizontal  layers  containing  constant  velocity 
gradients. 

If  is  the  velocity  gradient  between  depths  and  Zj^, 

ki.i  =  (Vi  -  Vi.i)/(Zi  -  Zi.i)  (50) 

Z=Zi_i  +  (V  -  Vi_j^)/ki_i  between  Zj^^  and  Zi.  (51) 

Combining  equations  (U9)  and  (5l)» 

Z  =  Zi_i  +  Vj.  (cos  ©  -  cos  Oi_i)/ki_i  between  and  Zi,  (52) 

Equation  (52)  clearly  represents  the  arc  of  a  circle  with  radius 
^r^^i-1*  follows  that  the  horizontal  component  of  the  circular  arc 
between  depths  Zi^^  and  Z^  is 

"3s^_l  =  Vj.  (sin  ei_^  -  sin  ©i)/ki_i  (53) 

where  Vj.  is  constant  between  successive  bottom  reflections  and  ki^^^  is 
constant  between  depths  Zi_2  and  Zi. 

Since  each  ocean  depth,  Zi,  has  only  one  corresponding  velocity, 

Vi,  a  ray  will  cross  Zi  with  identical  entrant  and  emergent  angles  of 
inclination,  ©i.  Reflection  from  the  level  ocean  surface  Will  result  in 
negative  entrant  and  positive  emergent  angles  of  inclination  which  are 
of  equal  magnitude.  Therefore,  substituting  equation  (53)  into  (U7), 
the  horizontal  distance  from  the  mth  bottom  reflection  to  the  nth  point 
on  a  ray  is 
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®i-l  - 

»  , 

where:  1.  All  velocity  gradient  layer  boundary  crossings  are  represented, 
2,  Surface  reflection  entrant  and  emergent  angles  of  inclination  are  of 
equal  magnitude  and  opposite  algebraic  sign.  3.  The  mth  bottom  re¬ 
flection  emergent  angle  of  inclination  is  linked  to  earlier  terminology 
by  931  =  ^i-1  obtained  from  equation  (50).  5.  (the 

ray ' s^'reversal  velocity  after  the  mth  bottom  reflection)  and  9^  are 
obtained  from  (the  velocity  at  reflection)  and  by  applying 
equation  (U9)  as  follows: 

V  ^i-1  -  ■  ^Bm  (55) 

^m  ■  cos  cos  9j^  ~  cos  9^jjj 

All  terms  in  expression  (4U)  are  now  available  except  the  stated 
partial  derivatives  of  and  0j^.  Combining  expressions  (5^)  and  (55) » 


n 

=  ^,^~*^(Vi-l  tan  9^_2^  -  tan 


Since  the  velocity  gradient  may  change  abruptly  at  a  layer  boundary  and 
since  the  ocean  surface  and  the  layer  boundaries  are  all  horizontal, 
discontinuities  in  the  calculation  of  9ASjj/3  9q  are  avoided  by  choosing 
3Zi/59  =  0,  where  is  the  depth  of  the  ith  point  on  a  ray.  Under 

this  condition,  equation  (56)  can  be  rewritten: 


aAS„  =  ^  (  V^•_-|  tan  9i_i  -  Vi  tan  9i  ^ 

t  - Jin  J  • 


39 

Since  velocity  is  constant  at  each  depth. 


3Vi_i/a9o  =  0,  3Vi/a9o  =  0, 
when  aZi_i/39o  =  0,  aZi/a9o  =  0, 
Employing  expressions  (50)  and  (58)* 


©o  =  0 


=  0 


(57) 


(58) 


(59) 
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Therefore,  from  expressions  (55)*  (57) »  (58),  and  (59), 

1  /  Vi_i  9  Oi_i  Vi  gQi  \ 

cos^  ©i  1  900  cos2  ©•  9©-/ 

^  ^  1  o/ 


1  /  1  3Qi.i _ 1  go. 

-^TTe’  i  ®i-l 


From  expression  (55), 

cos  ©i  _  Vi  cos  e^ni 

From  expressions  (58)  and  (61), 

-  sin  ©i  gQj  _  -Vj  sin  ©gni 

^  ®o  ^Bm  ^  ®0 


aSg  tan  ©i  a©Q 
Similarly, 

^  ^i-1  _  ®fim 

©i_i  a©o 

Substituting  equations  (62)  and  (63)  into  (60), 


®Bm  1  /  L. _ 1  \ 

99*,  ^sin  ©i_i  sin  ©j/ 

,  r 

Assuming  that  the  summation  increment  in  equation  (5^)  is  available,  it 
is  more  convenient  to  rewrite  the  above  as 


a^Sn  _  -Vj.  tan  0^^ 


Bm  Qj-l  -  sin  ©^ 

0  zCi  ^i-i  sin  ©i.i  sin  ©i 


23 


NAVWEP5  REPORT  T79S 


Deriving  3^S^/ in  a  manner  similar  to  equations  (56)  through  (64), 


■■  I  ■  -  •‘m  —  ■  r— 

3^0  <5  00 


sin  0j:_i  -  sin  0^^ 


■J0^/£_Ui_i  sin  Qi_i  sin  Oi 


The  above  equations  are  suitable  for  application  in  expression  (44).  It 
now  remains  to  determine  the  partial  derivatives  of  and  0jfi  in 
expressions  (44),  (64),  and  (65). 

Recall  that  equations  (39)  and  (4o)  display  a  relationship  between 
the  ray's  mth  bottom  reflection  entrant  and  emergent  ray  directions, 

(®Bra»  0K-l)  0m^ »  ."^^ere  0  is  the  ray's  angle  of  inclination 

and  0  is  the  angle  between  the  vertical  XZ  plane  and  the  vertical  plane 
containing  the  ray  path.  From  equation  (39), 

cos  04m  =  cos  0Bra  ^^Bm  _  pC  ^  sin  ©3^  ^^Bm  +  2  cos  03  ^^m3 


^«Bm 


^Bm 
cos  04nj 


to 

to  L 


’  2Cm3  tan  03^)  303^ 

5 


+  2 

J 


where,  from  equation  (4l), 


•Cm3  =  K-1  '^m-l  ^^n  1^^) 


tan^  +  tan^  ^ 


^  S>0Bin 

aSo”  03“  3  0o 


Deriving  from  equation  (4o), 


gV] 

'OS  ?9o 


-  2  tan  (T~r, 


cos2  0  dQ 


'^m-l  -  2Cj^3  tan 


.  sin  0„  T  -  2C_-5  tan  (7“„  (  ■  n  a0  t 

+  m-1 _ m3 _ m  /sin  0in-l  m-1  + 

(cos  0,^_i  -  2Cj„3  tan  9  ©o 


2  tan  ^  m3 


Employing  equation  (38)  in  the  above. 
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^0, 


cos  d 


Bm 


. V  ~  "  m  — 


cos*^  0m  3®o  ®Am 

0j^  (  sin  0j^_i  +  2  tan  i|<  U 

^  500  3«o 


a©. 


+  tan 


50jj^  cos  ©R„  cos  0, 


Bn  v>m 


cos  Q-i 


Bm 


“ 

(si 


^7  ^ 

in  0j„_2,  "*■ 


a©. 


+  2(tan  0n,  tan  'I' „  -  tan  <r-^)  1  . 

^Qo  J 


(68) 


Equations  (66)  through  (68)  describe  ^Qi^/^QQ  and  a0jn/^0Q  in 
terms  of  3  0Bm/®®o  ^0in-l/^®o«  However,  it  is  desirable  to  write 

a0Bm/^®o  terms  of  aQBj„.l/3©o  so  that  such  an  expression  can  be 
used  with  equations  (66)  through  (68)  to  establish  recursive  expressions 

for  9e6m/^®o  terms  of  9©]^m-l/ ^®o  90m-l'^  ^®o* 

The  desired  expression  for  9©Bm/9®o  obtained  from  the  following 

application  of  equation  (it8): 


cos  ©bj^  _  ^Bm  ®’Bm-l 
'^Bm-1 

where  0^,n_i  is  the  ray's  m-lst  bottom  reflection  emergent  angle  of  incli¬ 
nation  and  ©3jjj  is  the  ray's  mth  bottom  reflection  entrant  angle  of  incli¬ 
nation.  From  exijression  (58)  *3ind  the  above, 


-sin  ©Bm  _  ^3m  ^in  ©4n_i  9©^in.i 

*©0  ■  Vb^-I  ^®o 


^^Bm  -  •tan  egm-l  ^^3m-l  . 

a4i  tan  ©gjn  3  ©^  ( 69 ) 

Substituting  equation  (69)  into  (67), 

=  CmU  ^^m-1  *  Cm5  ^®5tn-l  (^0) 

a©o  500  500 
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where 


^raU  = 


cos  tan  cr^  -  sin  0j„_i  tan 


sin  01 


-  tan  0^^,;^ _ 

cos  0T3n.  (tan2  „  +  tan^  <7^  +  1) 


Noting  a  similar  derivation  of  and  5  0j^/9  0q,  and  substituting 

equations  (69)  and  (TO)  into  (66)  and  (68),  the  recursive  expressions 
for  the  partial  derivatives  of  0^jj^  and  0^^  are: 


(D 

II 

cos  0Bpi 

aeSm-l  +  2 

/cmH  30ni-l  +  CuiB 

90^1 1] 

300 

cos  06jn 

L  3^0 

V  aSo 

3  00  /J 

=  cos  Osn,  cos  0„,  ^0m-l  +  20^8  f Cmk  +  0^5  ^Q^m-1  \ 

3  00  cos  L  ’  3®o  \  ^«o  /- 

^^^m  ^  ^Bm  Fcjne  ^^Bm-1  ^  2  Ur^k  ^^m-1  +  Cj„5  ^ ^Bm-lj  (7I) 

3  00  cos  O^j^L  B0q  \  3  00  a  0c  / 

^  cos  0B;^  cos  0m  f a0m-l  +  2Cm8  ( ^mli  ^^m-1  +  Cm5  ^^^-1^1 

a0o  %m  L  900  V  P0O  500  /-J 

where 
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smd  from  equation 

tan^  ij<  m  +  tan^  a-“^  +  1 
GENERAL  RAY  TRACING  PROCEDURE 

The  spreading  loss  term  at  each  desired  point  on  a  ray  is  obtained 
from  expressions  (^2),  (5^)»  (55),  (64),  and  (65)  where,  before 

the  first  bottom  reflection, 

m  =  0  ^  ° 

Bm  =  0  axw  a©o  =  0  aY^^/  ao^  =  o  ax^/  =  o 

3YW  3!^o  =  0  efim  =  80  ^©Bm/  ^®o  =  ^  ^ 

=  ^o/cos  iSq 

When  the  ray  reaches  the  mth  bottom  reflection  point,  =  O^m* 

3V®«o  =  3XB„/ae„  aY„/a(io  -  sYj^/ao^ 

ax„/3«„  .  3XB„/a0o  aY„/a0o  =  aYBn/a0„ 

While  at  the  mth  bottom  reflection,  expressions  (44),  (54),  (64),  and 
(65)  are  up-dated  by  finding: 

1.  9]^m  ^m  equations  (39)  and  (40). 

2.  Vj-m  =  ^Bm  expression  (55). 

3.  aeV^®o.  ^«o.  ^8W^fi^o»  and  a0jj,/a0Q  from  expression 

(71). 

4.  ^®o»  ^^6m/^8o*  and  expres¬ 

sion  (28). 

After  the  mth  bottom  reflection,  spreading  loss  terms  at  desired  points 
on  the  ray  are  obtained  from  equation  (42)  and  up-dated  expressions 
(44),  (54),  (64),  and  (65). 

SPECIAL  EXAMPLES 

There  are  t..o  cases  of  particular  interest  which  greatly  simplify 
the  preceding  ray  tracing  procedure  for  calculating  the  spreading  loss 
term. 
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1.  The  initial  ray  path  lies  on  the  vertical  XZ  plane,  and  the  ocean 
bottom  is  described  such  that  it  is  orthogonal  to  the  XZ  plane. 

2,  The  ocean  bottom  is  a  horizontal  plane. 

CASE  I 


The  conditions  of  the  first  case  are  0©  =  0  and  O'-  0.  Substituting 
0  into  equation  (38a), 


sin  0jn  =  sin  ^ 

cos 

When  00  “  equation  (72)  yields 

K  =  '^m-l  =  ...  =  01  =  0o  =  0- 
Also  from  equation  (72), 


(72) 


(73) 


a0jji  cos  9  30 

cos  0jn  _  =  _ cos  0ni_i  +  sin  0^^,  ^  3  /cos  93^,  j  ^ 

0'^  3  0^  ■  a0^Vcos9^j,; 

From  expression  (73)  and  the  above. 


^0 


00=0 


_  cos  ©3^  a0j„.i 


Expressions  (72)  and  (73)  can  be  used  to  similarly  derive 


ilk) 


I  =  cos  ©Qjn  30ii,_i 

I  <»o='>  ■=<>=  «Bm  a»o 


Since  3  0pj/  SOq  ~  0  when  m  =  0,  the  above  equation  indicates  that  for 
any  m. 


a© 


m 


o 


0. 


(75) 


Since  3  0^/3  0^  =  1  when  m  =  0,  equation  (7^)  can  be  written 
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^  Ogja  cos  cos  Qgi  ^ 

00=0  ®Bin-l  ”* 

Restating  equation  (55), 


*  iu  - — 

cos 

where  Vj.  is  the  reversal  velocity  after  the  mth  bottom  reflection, 

Vgj^  is  the  velocity  at  the  reflection,  and  O^jj^  is  the  ray's  mth  bottom 
reflection  emergent  angle  of  inclination.  Equation  (55)  was  obtained 
from  equation  (^9),  from  which  a  similar  and  further  application  yields 


'rm_i  = 


’'^Bm-1 

®Bm-l 


''Bm 
cos  ©t 


where  ©gj^^  is  the  ray's  mth  bottom  reflection  entrant  angle  of  inclination. 
Combining  equations  (77)  and  (78), 

cos.©Bm^  Vr^  .  (79) 

cos  ©im 

Substituting  the  above  into  equation  (76), 

I  ^^m  ^^m-1 


Vr  V_  , 
^m  -^m-l 

^rm-1  Vrn_2 


Substituting  expressions  (73),  (75),  and  (80)  into  (^^), 


axn 

=  axSm 

o 

00=0  ° 

3Yn 

900 

00=0  ^®o 

^  ^^Bm  ^  gASn 

00=0 

^^n  _  ^  ^5m  ^^m  ^  “^n 

0  =0  ^  *^0 
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where  the  primed  partial  derivatives  are  found  from  ®^Bra/ 

a©o,  300,  and  3Y3m/30o,  using  expression  (28). 


Expression  (28)  is  an  ocean  bottom  boundary  condition  approximation 
which  can  be  simplified  by  employing  expression  (73)  and  ^=0. 


^^Bm 


tan  0 


Bm 


a©o 

00=0 

tan  0^^ 

_  ^^Bm 

300 

00=0 

3  0o 

S^Bm 

_  tan  ©gjg 

a«*o 

00=0 

tan  0^^ 

3Y^m 

= 

00=0 

WT 

ni- 

.1 

the  horizi 

tan  0^jjj  -  tan  \li 


m 


tan  03jjj  -  tan 


3  0^ 


tan  ~  m  \  9X 

tan  03j„  -  tan  ip 


■Bm 


(82) 


bottom  reflections.  Since  the  respective  coordinates  of  the  reflection 

points  are  (Xgj^.i,  ^Bm-l^  (%m»  ^Bm»  ^Bm)  ♦  ^3  (page  19) 

indicates  that: 


^Bm  -  ^Bm-l  +  -^^m-l  ^m-1 

^Bm  =  YBm-1  +  ^  ^m-1  (^m-1 

The  boundary  condition  discussion  preceding  expression  (28)  used  un¬ 
primed  and  primed  partial  derivatives  to  respectively  denote  before- 
and  after-reflection  rates  of  change.  Bearing  this  in  mind,  the  above 
equations  are  used  to  obtain 


^Xfim  = 

sin 

!^m-l 

3  0m-l 

+ 

cos 

0m-l 

3  A  S]^_i 

3  ©o 

3  00 

3  00 

3  00 

3YBm  - 

S'Yfim-l  + 

^  Sni_i 

cos 

0m-l 

3  0m-l 

+ 

sin 

0m-l 

<U 

(D 

o 

3  00 

3  00 

3  0o 

+ 

3  00 

^  Sm— 1 

cos 

^^m-1 

30m-l 
3  00 

+ 

sin 

0m-l 

Sin_i 

3  00 
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From  expressions  (73),  (75),  (80),  and  the  previous  equation. 


^^Bm 

3XiTn_i  ^ 

3AS„_i 

300 

00=0 

^00 

3  00 

3YBm 

(83) 

300 

00=0  ' 

9  00 

S^Bm 

^^Bra-l  4.  ^rni_i  ■^^m-1  • 

3  00 

o 

11 

o 

3  00  Vro 

Combining  expressions  (82)  and  (83), 


3XAm 

tan  ©Bn  / 

^tan 

04m  “  ^an  '1'  n  'j 

iZ-fatit 

300 

00=0 

tan  \ 

tan 

0Bm  -  tan  '*'n/ 

(~»o 

3Y4ni 

300 

o 

li 

O 

^00 

^^Bm 

3Y^m-l  4 

Vr™ 

^  171- 

o 

It 

o 

3  00 

< 

o 

Since  9X5ni/30o« 

3X;§in/  3  0, 

D,  and  3Y^j^/30o  = 

0  when  m  = 

expression  indicates  that  for  any  iri>0. 


3^0  / 


0,  the  above 


a  Or 


=  Om 


00=0 


+  d  ^  ^m-l) 

weo  / 


where 
and  Cjjji 


<9X]!^nj_l/ 900=0  when  m  =  1 
tan  Ogjjj  /  tan  0^jjj  -  tan 
tan  ^  tan  -  tan 


(8l») 


m 


a©. 


00=0 


0 
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m-/ 


From  expression  (8l)  and  above, 


9Y 

Te 


n 


o 


00=0 


=  0 


9Yn 

a0o 


AS 


J 


=  cos  0 


o 


AS. 


(85) 


(85a) 


Combining  equations  {85),  (85a),  and  (^2),  when  0o  ®  0  and  <J~'-  0,  the 
special  expression  for  the  spreading  loss  term  at  point  Pjj  on  a  ray  after 
m  bottom  reflections  is 


sin  0n  [  B  ©o 


A  S. 


where,  from  expression  (81), 


(86) 


(84). 


3^im  +  with  9X]^jjj/^0p  described  in  expression 

a^o  ^»o 
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Pertinent  terms  in  expressions  (86)  and  (81+)  are  defined  as  follows: 

is  the  velocity/  at  the  ray  source. 

Vj..  is  the  reversal  velocity  hetween  the  and  .j-i-lst  bottom 
'  reflections. 

AS.  is  the  horizontal  distance  between  the  and  j+lst  bottom 
'^reflections . 

ASjj  is  the  horizontal  distance  from  the  mth  bottom  reflection  to 
point  Pjj. 

is  the  mth  bottom  reflection  boundary  condition  correction. 

0Bj„  and  are  the  ray's  respective  mth  bottom  reflection  entrant 
and  emergent  angles  of  inclination. 

is  the  angle  the  ocean  bottom  tangent  plane  at  the  mth  reflec¬ 
tion  makes  with  the  horizontal  in  the  XZ  plane. 

It  is  interesting  to  note  that  Cjj^  =  p  if  the  boundary  condition  is  com¬ 
pletely  ignored,  as  was  done  in  reference  (^).  Under  this  condition, 
expression  (8lt)  yields 


which  is  substituted  into  expression  (86)  to  obtain 


Since  (SL)j^^  +  I^/In  from  equation  (l),  it  can  be  seen  that  the  above 
equation  agrees  with  equation  (l)  of  reference  (2). 

Expression  (86)  is  written  in  terms  of  ^Sj,  ASj^,  and  3ASjj/aOQ, 
where  ^  is  found  from  equation  (5^).  Similarly, 

BJ-h  I 

(sin  -  sin  ei)/ki_i  (8?) 
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3aSji/5©o  obtained  from  equation  (64) »  where  in  turn, 

is  described  in  expression  (71).  However,  when  {fg  “  ^  and  0,  a 
much  simpler  expression  for  3e^j^/a©o  can  be  found  by  first  substituting 
0~=  0  into  equation  (4l): 

Cj^3  =  (cos  tan  'I'm  -  tan  ©3^)  ’*'m 

Employing  expression  (73), 

Cm3  "I’m  -  ®Bm)  'I' m  ~  ^  * 

Substituting  the  above  into  equation  (39)*  where  ©gm  ®^m 

ray's  respective  mth  bottom  reflection  entrant  and  emergent  angles  of 

inclination, 

sin  =  sin  ©g^^  +  2  sin  cos  4'm  cos  ©Bm  -  2  cos^  il'm  sin  ©Bm 


=  sin  (2  if 

m 


®Bm 


) 


«Bm  =  ^ 


’*’m  ‘  ®Bm 


when  0Q  =  0 


(88) 


From  equations  (88)  and  (69), 


B)o=0 


^^Bm  _  Q^nui  ^Q^m-1 

3©o  tan  ©Bm  ^ 


Noting  that 


(89) 


®Bm  =  ®o  ^®Bm/  ^®o  =  1  m  =  0, 


(90) 


equation  (89)  indicates  that  for  any  m  >0, 


nn-/ 


II 

_ 1 

0 

II 

0 

TS. 

t3.n, 

tan  ©Bj+3^ 


tan  ©Q 

©4m 


07 


tan  ©^ 
tan  ©B 


(91) 


Substituting  equations  (90)  and  (91)  into  (64), 
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tan 

tan  Ob 


Oj.l  -  sin  0^ 
sin  sin  Oj 


[92) 


where  the  continued  product  equals  one  when  m  =  0. 

When  0Q  =  0  and  <7^=  0,  the  spreading  loss  term  at  each  desired  point 
on  a  ray  after  m  bottom  reflections  is  obtained  from  expressions  (86), 
(5H),  and  (92).  Just  as  the  ray  reached  the  mth  reflection, 

0^  =  ©Bm  AS„  =  aAS^/ao^  = 

with  expressions  (86),  (5^),  and  (92)  up-dated  through  reflection  by 
finding: 


®^m  “  2  -  Og^  from  equation  (88). 

2.  =  Vb^/cos  O^jjj  from  equation  (77). 

3.  expression  (84), 


The  validity  of  spreading  loss  expression  (86)  is  partly  substantiated 
by  its  agreement  with  equation  3B-42  of  reference  (l)  when  no  bottom 
reflections  are  involved.  Under  this  condition,  m  =  0  and  9X^/9  0^  =  0, 
which  are  used  in  expression  (86)  to  obtain 


(SL)  = 

Vo 

/  a  bn] 

li 

sin  Ojj 

J  1 

cos  0Q 

/as^ 

• 

sin  ©n 

\ 

(93) 


Let  be  the  horizontal  distance  from  ray  source  to  point  Pj^.  Then 
before  the  first  bottom  reflection,  Sjj  =  AS^,  previously  defined  as  the 
horizontal  distance  from  the  mth  bottom  reflection  to  point  Pj^.  There¬ 
fore,  from  equations  (l)  and  (93),  the  spreading  loss  term  at  point  P^ 
on  a  ray  before  any  bottom  reflections  is 


(sl)-3  =  ^ 


sin  0|^  Sj^  9  Sjj 
cos  Oq  aOo 


which  agrees  with  equation  3B-42  of  reference  (l). 


(94) 
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CASE  II 

The  ocean  bottom  is  a  horizontal  plane  when  its  tangent  plane  in¬ 
clination  angles,  and  are  zero.  A  ray  path  leaves  its  source 
with  direction  (Qq,  0o^»  'where  is  the  ray's  initial  angle  of  inclina¬ 
tion  and  0o  is  'the  angle  between  the  vertical  XZ  plane  and  the  vertical 
plane  containing  the  ray  path  before  any  bottom  reflection.  When  =  0 
and  (T^-  0,  the  environment  is  everywhere  symmetric  about  the  vertical  Z 
axis  passing  through  the  ray  source.  Therefore,  the  family  of  ray  paths 
with  a  common  source  and  a  common  initial  angle  of  inclination  will 
have  identical  physical  characteristics  at  a  given  horizontal  range,  re¬ 
gardless  of  the  value  of  0q.  The  simplest  representative  ray  to  study 
is  that  in  which  0q  =  0. 

Since  the  conditions  of  Case  I  were  cr'=  0  and  0o  =  0  and  the  condi¬ 
tions  of  this  horizontal  bottom  case  are  '/'  =  0,  =  0,  and  chosen  0^  =  0, 

Case  I  spreading  loss  expression  (86)  can  be  used  to  derive  the  spreading 
loss  expression  for  this  case.  Setting  ^■=  0  in  expression  (81|), 


34SJ 


00=0  ^ 

I 

which  is  substituted  into  expression  (86)  to  obtain 


'o  00=0 


^  3ASn 
,3  ©o  3  ®0 


If  Sjj  is  the  horizontal  distance  from  ray  source  to  point  Pj^, 


.1 


It  has  been  consistently  assumed  that  3Zi/ dQo  =  0,  where  Z^  is  the 
depth  at  any  given  point  on  a  ray.  Therefore,  since  reflection 
points  on  the  horizontal  ocean  bottom  satisfy  the  above  condition, 
equation  (97)  can  be  rewritten 
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in~l 

J-  o 

When  the  ocean  hottom  is  horizontal, 
simplified  by  noting  the  equality  of 


(97a) 


spreading  loss  expression  (86)  is 
equation  (95)  and  the  above. 


ax. 


ae. 


00=0 


_ffn 


(98) 


A  horizontal  bottom  also  simplifies  the  relationship  between  entrant 
and  emergent  inclination  angles  of  a  ray  at  reflection.  Setting  ^  =  0 
in  equation  (88), 


“  “  ®Bm 


(99) 


Employing  the  above  in  equation  (79), 
Vr„/V-  =1 


It  follows  that 


“  ^rm-l  “  •••  " 

Equations  (96)  and  (lOO)  indicate  that 


(100) 


/rt-'/ 


AS: 


+  Vr 

^m  a 


Sn  = 


Vq  Sn 

cos  0Q 


Substituting  equation  (98)  and  the  above  into  expression  (86),  when  the 
ocean  bottom  is  a  horizontal  plane,  the  special  expression  for  the 
spreading  loss  term  at  point  Pn  anywhere  on  a  ray  is 
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which  happens  to  agree  with  equation  (9^). 


(101) 


An  evaluation  of  is  obtained  from  equations  (96),  (87),  and  (5^+). 
=  X  ^rj  /  ^ (sin  -  sin  ^  (sin  -  sin 

J~o  \i:  f+sj  J  t= 


Employing  equations  (99)  and  (lOO)  in  the  above, 

n 

Sn  =  VrQ2^(sin  Qi.i  -  sin  ei)/ki_i  (102) 

»■*/ 

where  entrant  and  emergent  angles  of  inclination  at  reflections  are  of 
equal  magnitude  and  opposite  sign.  A  similar  handling  of  equations 
(97a),  (92),  (99) »  and  (lOO)  produces 


=  “  ^ro  "tan  ©q  sin  9^,1  -  sin  Qj  j  (103) 

300  ^i-1  ®i-l  ®i/ 

/■=/ 

Equations  (lOl),  (102),  and  (103)  provide  a  detailed  expression  for  the 
spreading  loss  term  at  the  nth  point  on  a  ray  in  an  environment  con¬ 
taining  a  horizontal  ocean  bottom. 

CONCLUSION 

The  general  ray  tracing  procedure  on  page  27,  the  special  ray  trac¬ 
ing  procedure  on  page  35»  and  equations  (lOl),  (102),  and  (103)  provide 
spreading  loss  descriptions  which  are  developed  for  use  with  different 
ocean  bottom  requirements.  Each  approach  is  designed  for  efficient 
application  on  a  high-speed  digital  computer. 
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Appendix 

SPECIAL  APPLICATION  OF  THE  AREA  OF  A  QUADRILATERAL 


THEOREM 

The  area  of  a  quadrilateral  is  equal  to  eight  times  the  area  of  any 
triangle  whose  vertices  are  the  midpoints  of  two  ad.'jacent  sides  of  the 
quadrilateral  and  the  point  of  intersection  of  the  two  lines  Joining 
opposite  midpoints. 

PROOF 

Construct  quadrilateral  ABCD  with  diagonal  d  between  points  B  and  D, 
and  with  side  midpoints  E,  F,  G,  and  H  as  shown  in  Fig.  lA. 


Drop  perpendiculars  h  and  i  from  respective  points  A  and  C  to  diagonal 
d.  From  Fig.  lA, 

Area  of  quadrilateral  ABCD  =  hd/2  +  id/2  =  (h  +  i)  d/2.  (l) 

Construct  lines  joining  side  midpoints  E,  F,  G,  and  H  as  shown  in 
Fig.  2A.  From  triangle  similarities  (examples:  GCH  and  BCD,  EAF  and 
DAB) , 

EFGII  is  a  parallelogram  with  base  d/2  and  altitude  (h/2  +  i/2).  (2) 

Construct  the  diagonals  of  parallelogram  EFGK  with  their  intersection 
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at  P. 

Theorem;  The  diagonals  of  a  parallelogram  bisect  each  other.  (3) 

From  expressions  (2)  and  {3),»  triangles  FPE  and  HPG  are  congruent,  with 
base  d/2  and  altitude  (h  +  i)/U.  Therefore, 

Area  of  triangle  FPE  =  (h  +  i)d/l6.  {h) 

Combining  equations  (l)  and  (4), 

Area  of  quadrilateral  ABCD  equals  8  times  area  of  triangle  FPE.  (5) 

By  dropping  a  perpendicular  from  point  E  to  diagonal  FH,  it  is  seen  that 
the  arWs  of  triangles  FPE  and  EPH  are  equal,  since  they  have  equal 
bases  and  a  common  altitude.  It  can  be  similarly  shown  that 

Triangles  FPE,  EPH,  HPG,  and  GPF  all  have  equal  areas.  (6)^ 

Expressions  (5)  and  (6)  complete  the  proof. 
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